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We address the question of the dependence of the bulk free energy on boundary conditions 
for the six vertex model. Here we compare the bulk free energy for periodic and domain wall 
boundary conditions. Using a determinant representation for the partition function with 
domain wall boundary conditions, we derive Toda differential equations and solve them 
asymptotically in order to extract the bulk free energy. We find that it is different and bears 
no simple relation with the free energy for periodic boundary conditions. The six vertex 
model with domain wall boundary conditions is closely related to algebraic combinatorics 
(alternating sign matrices). This implies new results for the weighted counting for large 
size alternating sign matrices. Finally we comment on the interpretation of our results, in 
particular in connection with domino tilings (dimers on a square lattice). 
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1. Introduction 

The six vertex model is an important model of classical statistical mechanics in two 
dimensions. The prototypical model is the ice model, which was solved by Lieb [1] in 1967 
by means of Bethe Ansatz, followed by several generalizations [2], [3], [4]. The solution 
of the most general six vertex model was given by Sutherland [5] in 1967. The bulk free 
energy was calculated in these papers for periodic boundary conditions (PBC). A detailed 
classification of the phases of the model can be found for example in the book [6] (see also 
the more recent work [7] on anti-periodic boundary conditions). 

Earlier, in 1961 Kasteleyn, while studying dimer arrangements on a quadratic lattice, 
expressed doubts on the independence of the bulk free energy on boundary conditions [8] . 
For more on dimer arrangements, see [8], [9] and [10]. Interest in this subject was renewed 
with recent work on domino tilings (which are equivalent to dimers on a square lattice) of 
an Aztec diamond [11], [12], demonstrating a strong effect of the boundary on a typical 
domino configuration (see also [13]). Dimers (or domino tilings) can be considered as a 
particular case of the six vertex model, and therefore a natural question is to investigate 
the effect of boundary conditions on the thermodynamic limit of the six vertex model. 

Independently of this, new boundary conditions of the six-vertex model, the so-called 
domain wall boundary conditions (DWBC), were first introduced in 1982 [14] (we shall 
define them in detail below). An important recursion relation for the partition function 
was discovered in this paper. Later these recursion relations helped to find a determinant 
representation for the partition function of the six vertex model with DWBC [15], [16]. 
The determinant representation simplifies somewhat in the homogeneous case. In this 
case the partition function satisfy Toda differential equation [17]. In this paper we use this 
differential equation in order to calculate the bulk free energy for domain wall boundary 
conditions. 

Let us all mention that there is a one to one correspondence between arrow config- 
urations in the six vertex model with DWBC and Alternating Sign Matrices (ASM) [18]. 
This mapping was used in order to count the number of ASM. More on ASM can be found 
in [19] and [20]. 

The plan of the paper is as follows. In Section 2 we define the six- vertex model with 
domain wall boundary conditions, and derive the determinant representation for the parti- 
tion function. In Section 3 we derive Toda differential equation for the partition function. 
In Section 4 we consider the thermodynamic limit; we derive the explicit expression of 
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the bulk free energy in the ferroelectric and disordered phases, and compare it with PBC. 
Finally, in Section 5 we comment on the connection of our results with other subjects 
(ASM, domino tilings, height model) and conclude this discussion in Section 6. 



2. Determinant representation of the partition function of the six-vertex model 

In this section we shall define the inhomogeneous six-vertex model with domain wall 
boundary conditions, and rewrite its partition function as a determinant. We will then 
particularize our formula to the homogeneous case. 
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Fig. 1: A configuration of the inhomogeneous six- vertex model with domain 
wall boundary conditions. 



First we define the configurations of the model. They are given by assigning arrows to 
each edge of a N x N square lattice (see Fig. 1). The "domain wall" boundary conditions 
correspond to fixing the horizontal external arrows to be outgoing and the vertical external 
arrows to be incoming. The partition function is then obtained by summing over all possible 
configurations: 

TV 

z= Yi n ^^fe (2.1) 

arrow configurations i,k=i 

where the statistical weights Wik are assigned to each vertex of the lattice. Since we are 
considering an inhomogeneous model, we need two sets of spectral parameters {Aj} and 
{/i/s} which are associated with the horizontal and vertical lines. The weight Wik depends 
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on the arrow configuration around the vertex {i,k) and is given by 



a{Xi,iJ,k) 



->- 



Wik = { b{Xi,iJ,k) 



->- 



->- 



(2.2) 



->- 



->- 



(all other weights are zero) where the functions a, b, c are chosen as follows: 

a(A, fi) = sinh(A — — 7) 
6(A, //) = sinh(A — + 7) 

c(A, /i) = sinh(27) 



(2.3) 



Here 7 is an anisotropy parameter which does not depend on the lattice site. The partition 
function is therefore a function of the 2A'" spectral parameters and we shall denote it by 
^iv({Aj},{//fc})- 

The model thus defined satisfies the following essential property (Yang-Baxter equa- 
tion) shown on Fig. 2. The vertex with diagonal edges is assigned weights (the so-called 
R matrix) which are the same as the usual weights, up to a shift of the difference of the 
spectral parameter. Here we shall not need the explicit expression of the R matrix. 




Fig. 2: Yang-Baxter equation. Summation over arrows of the internal edges 
is implied, whereas external arrows are fixed. 



We shall now list the following four properties which determine entirely Ziv({Aj}, {/ife}) 
and sketch their proof (for a detailed algebraic proof the reader is referred to [16]): 
a) Zi = sinh(27). 

This is by definition. 
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b) -Z'jv({Ai}, {Hk}) is a symmetric function of the {A^} and of the 

It is sufficient to prove that exchange of /i^ and (for any i) leaves the partition 
function unchanged. This can be obtained by repeated use of the Yang-Baxter property: 
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(2.4) 



where R^^ = Rn is the appropriate entry of the R matrix; and similarly for the {A^}. 

c) ZAr({Ai}, {fik}) — e~*^^~^^''^*Pjv-i(e^^') where Pn-i is a polynomial of degree N — 1, 
and similarly for the /x^. 

Let us choose one configuration. Then the only weights which depend on Aj are the N 
weights on row i. Since the outgoing arrows are in opposite directions, at least one of the 
weights must be c. Therefore there are at most N — 1 weights a and b, and the product 
of all weights is of the form e~'^^~^^^*PAr-i(e^^*). This property remains of course valid 
when we sum over all configurations. 

d) Zn{{K}i {A*fe}) obeys the following recursion relation: 

^iv({Az},{yUfc})|A,-Mi=7 =sinh(27) JJ sinh(Aj - + 7) 



l<fc<Ar 
k^l 



JJ sinh(Aj - /// + 7) ZN-i{{Xi}i^j, {lJ.k}k^i) 



(2.5) 



l<i<N 



Because of property b), we can assume that j = I = 1. Since \k — l^-l = 1 implies 
a{\j — Hi) = 0, by inspection all configurations with non-zero weights are of the form 
shown on Fig. 3. This immediately proves Eq. (2.5). 
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Fig. 3: Graphical proof of the recursion relation. 



It is easy to see that the four properties a), b), c) and d) characterize entirely 
■^7v({Aj}, {//fc}). This is enough to prove that Zjv({Ai}, {//fc}) has the following deter- 
minant representation [16]: 



^7v({Aj},{//fc}) 



ni<i,fc<7v sinh(Ai - /Xfc + 7) sinh(At - jik-j) 
ni<i<,<Ar sinh(Ai - \j) ni<fe<z<iv sinh(//fe - 



det 

l<i,k<N 



sinh(27) 



sinh(Ai - iik + l) sinh(Ai - /ik - l) 



(2.6) 



Indeed, one can check that this expression satisfies the four properties listed above. 

The expression (2.6) might seem singular when two spectral parameters Aj and Xj 
coincide (and similarly for the /ik)', but in fact the pole created by the factor sinh(Aj — A^) 
is compensated by the zero of the determinant due to the fact that two rows are identical. 
Therefore, particular care must be taken when considering the homogeneous limit where 
all the Xi are equal (and all the fi^)- This limit was studied in detail in [16], and we shall 
simply summarize the result of the calculation. Let us call t the common value of A^ — 
for all i and k. When the Xi are close to one another one must Taylor expand the function 



0(t)^ 



sinh(27) 



sinh(t + 7) sinh(t — 7) 
which appears in the determinant. This leads to the following expression: 



(2.7) 



(2.8) 



3. Determinant representation and Toda chain hierarchy 

We shall now investigate the properties of the determinant which appears in Eq. (2.8), 
and for which we introduce the notation 



(3.1) 



with 



d" 



(3.2) 



Let us write down the bilinear Hirota equation satisfied by the r^v- For completeness, 
we recall that they are a consequence of Jacobi's determinant identity: 



(3.3) 



The large squares represent a given matrix, and the shaded regions are the sub-matrices 
whose determinants one must consider. Applying it to tn+i (up to a re-shuffling of the 
rows and columns), we find [17]: 





TnTn - = Tjv+lTjV-1 V7V > 1 



(3.4) 



where primes denote differentiation with respect to t. This is supplemented by the initial 
data: tq = 1 and ti = 0. Equivalently, we have: 



(logTjv) = 2 ViV > 1 



(3.5) 



'N 



which is the form of the equation that we shall use. 

Note that if we introduce the combinations e"^^ = tn/tn-i, N >1, Eq. (3.5) implies 
for the (fiN' 

ip'^ = QfN+i-fN _ qVn-'pn-i V7V > 2 (3.6) 

and (fi = _ These are the usual Toda (semi- infinite) chain equations [21], [22], [23]. 

Another possible form is 



V'^(r = -^C'Mive^- N>1 



(3.7) 



M 



with t/jjv = v'A'^+i ~ fN and Cmn (M, N > 1) the Cartan matrix of the semi-infinite 
diagram Aoq. 

This suggests a connection with the Toda chain hierarchy [24], [25], [26], [27], [28], 
[29]. Indeed, let us mention that given a Hankel matrix {mi+k-2) ^ i-e. whose entries only 
depend on i + k - the can be made to depend on a set of parameters {tq}q>i in such a 
way that the determinants r^r become r-functions of the whole Toda (semi-infinite) chain 
hierarchy [25]. Namely, one must choose 

mn{{tq}) ^ J dp{x) x"" e^''>^*'''^ (3.8) 

where dp{x) is an arbitrary measure"*^ (in the matrix model context [25], the tq are the 
coefficients of the polynomial potential). Here, we are in the simplest situation where only 
one parameter ti = t is allowed to evolve. We immediately check that Eq. (3.8) implies 
that rrinit) = -^Tnoit), which is consistent with Eq. (3.2). 



4. The thermodynamic Hmit 

We shall now consider the thermodynamic (i.e. large N) limit of the expression (2.8) 
in the various regimes of the six-vertex model. For that we shall use the Hirota equation 
in its form (3.5). 

When A?" 00 it is expected that the partition function behaves in the following way: 

log ZN{t) = -N^F{t) + 0{N) (4.1) 

where F{t) is the bulk free energy (we shall always set the temperature ksT = 1). Our 
main goal is to compute explicitly F{t). 

Comparing the expected asymptotic (4.1) with the exact formula (2.8), we find that 
the determinant rjv must be of the form 

rj,= (]in\\ e^'^(*)+0(^) (4.2) 

where 

f{t) = -F{t) - log(sinh(t + 7) sinh(t - 7)) (4.3) 

^ This must be considered as a formal expression; e.g. the measure may not necessarily be 
smooth or positive. 
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We now want to substitute the expansion (4.2) into the equation (3.5). For that we need 
to assume that the sub-dominant corrections to the bulk free energy vary slowly as a 
function of N; we shall discuss the validity of this assumption below. We then find that 
the expansion is consistent since both left and right hand sides of (3.5) turn out to be of 
order N^. The resulting equation for / is: 

/" = e'f (4.4) 

This is an ordinary second order differential equation, which can be readily solved. The 
general solution depends on two parameters a and to: 

e^^*^ = — TT-^ ^ (4.5) 

smh(a(t-to)) 

If the weights are chosen to be real, then the free energy should be real and this implies 
that a must be real or purely imaginary. 

So far everything we have done was independent of the particular form of the function 
(f){t) and therefore independent of 7. In order to fix the two constants in (4.5), we must 
now discuss separately the different regimes of the six-vertex model. Let us recall that the 
latter are usually distinguished by the value of the parameter (cf Eq. (8.3.21) of [6]) 

+ b- - 

A = ^ , (4.6) 
The weights a, 6, c were defined in Eq. (2.3) (with A — = t). In this parameterization, 

A = cosh(27) (4.7) 

Jf..l. Ferro-electric phase: A > 1 

This corresponds to the parameters 7 and t real; we recall that the weights are given 

by 

a = sinh(t — 7) b = sinh(t -|- 7) c = sinh(27) (4-8) 

with I7I < t. This is the so-called ferroelectric phase. In the case of periodic boundary 
conditions, it is known that the system is frozen in its ground state configuration, in which 
all arrows are aligned: if a > 6 all arrows point up and to the right or down and to the 
left, whereas if 6 > a they point up and to the left or down and to the right. The domain 
wall boundary conditions do not allow all arrows to be aligned: the ground state will 
instead take the form of Fig. 4. However at leading order in the large N limit, this does 
not affect the free energy, and we expect to find the same result as for periodic boundary 
conditions. 
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Fig. 4: Ground state configuration of tfie ferroelectric phase (for 6 > a; the 
case o > 6 is obtained by taking the mirror image). 

Indeed, it is easy to see that the relevant solution of Eq. (4.4) is 

1 



sinh(t — I7I) 

and therefore the bulk free energy takes the form 

g--F'(t) _ sinh(t + I7I) = max(a, 6) 
in agreement with the case of periodic boundary conditions. 

Jt..2. Disordered phase: — 1 < A < 1 

In this regime, it is customary to make the following redefinitions: 

7' = z(7 + i7r/2) 

t' = i{t + m/2) 
and divide all the weights by z, so that they take the form: 

a = sin(7 — t) b = sin(7 + t) c = sin(27) 



(4.9) 



(4.10) 



(4.11) 
(4.12) 

(4.13) 



and A = — cos(27). Using symmetry considerations, one can always assume < 7 < 7r/2. 
We only consider the region |t| < 7 (where the weights are positive). 

Taking into account these redefinitions, the partition function becomes: 



ZN{t) = 



(sin(t + 7) sin(t-7))^ 



(nS"!)' 



det 

l<i,fe<Ar 



]i+fc-2 



(4.14) 



with a redefined (f){t) = sin(27)/(sin(t — 7) sin(t + 7)); the determinant tn still satisfies 
Eq. (3.4) and f{t) defined by (4.3) is still a solution of Eq. (4.4). ^ 

Let us mention that the partition function has been computed exactly [30] at three 
particular values of the parameters: t = 0, 7 = 7r/6, /t/4 and n/3. In all three cases the 
expansion (4.1) and the assumption of smoothness of the sub-dominant corrections (which 
is necessary to derive the ordinary differential equation (4.4)) can be checked exactly. We 
have also checked it numerically for a variety of values of t and 7. 

We must now select the appropriate solution (of the form (4.5)) of the Eq. (4.4). Let 
us first assume that \t\ < 7 (this is the only physical region, i.e. where all the weights 
are positive). It is easy to check that f{t) must be an even function of t. The only even 
solution of Eq. (4.4) is 

em = (4.15) 
cos (at) 

where a remains to be determined. Note that this implies for F{t) 

e-^(*) = sinh(7 - t) sinh(7 + t) — ^ (4. 16) 

cos (at) 

We must then use the boundary condition given by \t\ = ±7. At these values one can 
compute directly ZN{t). Indeed the only non-zero configurations are of the form of Fig. 4, 
and we find 

ZAr(t = ±7) = sin(27)^' (4.17) 

and therefore e~^^^^ = sin(27). Since the prefactor in (4.16) vanishes when \t\ = 7, we 
conclude that a must be chosen in such a way that cos(at) is non-zero for \t\ < 7, but 
vanishes as \t\ = 7. This uniquely determines a to be: a = We obtain the final 
expression 

e"^^*) = sin(7 - t) sin(7 + t) ^^^7 ^ (4-18) 
^' ' ^' ^cos(7rt/27) ^ ^ 

As a consistency check, one takes the limit t ±7 and finds e"''^*-*^ = sin(27), as it should 
be. Also, note that for 7 = 7r/4, where the partition function is known and independent 
of t, one finds indeed that e"^*^*-* = 1. 

For further checks, let us set t — 0; a more standard normalization of the weights is 

then 

a = 6=1 c = 2cos7 (4.19) 



^ Note that the sign is unchanged in Eq. (4.4); this is the combined effect of the "Wick rotation" 
of t {t —>■ it) and of dividing all the weights by i (e-^ —ie^). 
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and the bulk free energy becomes 



_p TT sm 7 



(4.20) 



At 7 = tt/6, it/ 4, 7r/3, the values predicted by (4.20) coincide with the large N limit of 
the expressions of [30]. Also, this fits perfectly with some numerical computations of the 
determinant we have performed. 

We can compute the bulk energy (energy per unit site), which turns out to be 



E = (cot 7 — 1/7) cot 7 log(2cos7) 



(4.21) 



Fig. 5 shows the comparison with Monte-Carlo simulations. The agreement is also very 
good. 
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-0.5 



0.5 



Fig. 5: Energy E as a. function of the anisotropy A. The curve is given by 
Eq. (4.21), whereas the diamonds are the results of Monte-Carlo simulations 
on lattices of size N = 64. 



Finally let us mention that there seems to be no simple relation between the PBC 
and DWBC bulk free energies: from an analytic point of view, the DWBC free energy is 
an elementary function, whereas the PBC free energy is given by a non-trivial integral. 
Furthermore, the DWBC free energy is always greater then the PBC free energy, even at 
infinite temperature (A = 1/2), see Fig. 6. 
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Fig. 6: Bulk free energies for PBC and DWBC as a function of A. 
4.3. Anti-ferroelectric phase: A < — 1 

In this phase, the smoothness assumption of the sub-dominant corrections to the bulk 
free energy is not satisfied, as can be clearly seen numerically. The ratio Zjv+i^Ar-i/-Z^|r 
does not converge in the large N limit but instead has a pseudo-periodic behavior rem- 
iniscent of the one-matrix model with several cuts [31], and slightly more sophisticated 
methods are needed to analyze the large N limit. We leave this to a future publication. 

4.4- Phase transition at A = 1 

If the Boltzmann weights depend on a parameter (e.g. temperature), it is known that 
with periodic boundary conditions, the system undergoes phase transitions as A crosses 
±1. Let us use the expressions of the bulk free energy found above to clarify what happens 
in the case of domain wall boundary conditions. 

Here we shall consider the transition from ferroelectric (low temperature) to disordered 
(high temperature) regime, that is from A > 1 to A < 1. The parameter that plays the 
role of deviation from criticality T — Tc can be defined as 



We assume that b > a (the case a > b can be treated similarly), and re-scale the weights 
so that 6=1. With this convention, we simply have in the ferroelectric phase: 



T - Tc = 1 - A 



(4.22) 



A > 1 



(4.23) 
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(cf Eq. (4.10)). Let us now consider A — > 1 . The weights are 

sin(7-t) sin(27) 



sin (7 + t) sin (7 + t) 



(4.24) 



with 7 = 7r/2 + e, t = 7r/2 + ex; x must be kept fixed as e — > 0. Note that A = cos(2e), so 
that 

T - Tc oc (4.25) 
Expanding the expression (4.18) for the free energy, we obtain: 

e-^ = 1 - ^i^-lil!e3 + o(,4) (4 26) 

OTT 

Comparing (4.23) and (4.26), we find a second order ph.a,se transition, with a singular part 
(T — Tc)3/^ corresponding to a critical exponent a = — 1/2. This is to be contrasted with 
the first order phase transition that occurs in the case of PBC. Let us however emphasize 
that the difference of orders of the phase transitions is not that significant, since the phase 
transition is somewhat special (in the case of PBC, the correlation length jumps from zero 
for A > 1 to infinity for A < 1). 



5. Some equivalences 

We shall now review some alternative interpretations of the partition function of the 
six-vertex model with domain wall boundary conditions; these equivalent formulations will 
shed some light on the property of dependance on boundary conditions that was found. 

5.1. Alternating sign matrices 

Six-vertex model arrow configurations with domain wall boundary conditions on a 
N X N lattice are in one-to-one correspondence with alternating sign matrices (ASM) of 
size N, that is square matrices with entries or ±1 such that each row and column has 
an alternating sequence of -|-1 and —1 (zeroes excluded) starting and ending with a -|-1. 
Recalling that there are 6 weights which we shall call ai, 02, bi, 62, ci, C2 in the order 
shown in Eq. (2.2), the correspondence goes as follows: given a six- vertex configuration, 
assign a to each vertex a or 6 and -|-1 (resp. —1) to each vertex ci (resp. C2). One can 
show that this map is bijective, and therefore, the number of ASM is exactly equal to the 
partition function considered before with a = b = c = 1. For our purposes, let us define a 
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refined counting of ASM (a;-enumeration in the language of [30] ) by assigning a weight x to 
each entry —1 of the ASM. The resulting quantity A{N, x) is still related to the six-vertex 
model; indeed, one can easily show that 



A{N, x) = x-^l'^Z^ia = b=l,c= y^) 



(5.1) 



If < a; < 4, one can set x = 4cos^7 and the weights are of the form (4.19). One can 
then prove (extending slightly the asymptotic expansion found in 4-2) that 



log A{N, x) = N"^ log 



TT sm7 



N 

— \ogx + 0{\ogN) Va;e[0,4] (5.2) 
2 



Though this equivalence does not directly provide any useful insight into the issue adressed 
in this paper, the result (5.2) itself might be of some mathematical interest. 

5.2. Tilings of the Aztec diamond 

A more illuminating equivalence is that of domino tilings (i.e. dimers on a square 
lattice in a dual description) and the six-vertex model at A = - both models are well- 
known to describe essentially one Dirac fermion. This is illustrated on Fig. 7. Since each 
vertex of type ci has 2 possible corresponding domino tiling configurations, one must assign 
it a Boltzmann weight of 2 in order to count each domino tiling exactly once; however with 
most boundary conditions there are as many vertices of type ci and C2, and therefore one 
can give them both a weight of -\/2 instead, which leads to the values a = 6 = 1, c = \/2 
of the parameters. 



flj - — >- 



->- 



Fig. 7 : Correspondence between vertices of the six- vertex model and small 
patches of a domino tiling. 
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a) b) 

Fig. 8: a) A configuration of the six-vertex model with DWBC, and b) one 
possible corresponding tiling of the Aztec diamond. 

The more precise statement is that the number of domino tilings of the Aztec diamond 
(see [11]) is equal (up to a small known prefactor) to the partition function of the six- vertex 
model with domain wall boundary conditions ata = 6 = l, c = \/2, see Fig. 8. 

Since this a local correspondence of configurations it entends to all correlation func- 
tions. Also, introducing some weights for the local tiling patterns amounts to changing 
the weights a, b, c, but always in such a way that A remains zero. 

These tilings have been an object of interest for mathematicians, see in particular [12], 
[18]. The "arctic circle theorem" [11] shows that as the size of the system grows large, the 
domino configurations become frozen outside the circle inscribed inside the diamond, and 
remain disordered but still heterogeneous [12] (i.e. non translationally invariant) inside the 
circle. These statements have a straightforward equivalent in the six-vertex language: we 
expect the one-point functions of the six- vertex model at A = with DWBC to be non- 
constant (following a similar pattern as the tilings), and presumably a similar behavior at 
A 7^ 0. This gives a qualitative understanding of the dependence of the bulk free energy 
on the boundary conditions. 

5.3. Height model 

The general eight-vertex is well-known to be equivalent to a class of height models 
(SOS/RSOS model). In the case of the six- vertex model, there is a particularly simple 
equivalence which goes as follows: given a six-vertex configuration, integers are assigned 
to each face of the lattice in such a way that going from one face to a neighboring face, 
the number is increased by one if the arrow in between goes right (and so, is decreased 
by one if it goes left). Conservation of the arrows ensures consistency of this procedure. 
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The Boltzman weights of the model are simply the weights of the original six- vertex model 
expressed in terms of the new height variables. 

This equivalence is particularly interesting because it gives us a simple intuitive ex- 
planation of the lack of thermodynamic limit due to boundary conditions; this was studied 
in detail and proven rigorously in the case of tilings (A = 0), see [13]. Let us consider 
the domain wall boundary conditions and translate these into the language of our height 
model. They are fixed boundary conditions for the heights, of the form: 

1 ••• N-1 N 

1 N-1 

N-1 1 
N N-1 ■■■ 1 

where we have fixed arbitrarily the upper left height to be zero. 

In the thermodynamic limit A?" — > oo, let us define the rescaled coordinates on the 
square lattice to he x = k a, y = i a where a = 1/N is the lattice spacing. The heights hij 
are supposed to renormalize, according to standard lore, to a free massless bosonic field: 

hij (j){x, y) (5.3) 

However, it is reasonable to assume that in order to have a proper thermodynamic limit, 
the boundary conditions must be well-defined in terms of the limiting field cj). In the case 
of DWBC, one finds that the boundary conditions become (t){x, 0) = x/a etc, which do not 
have a limit as a ^ 0; in particular, the variations of (j) on the boundary diverge. More 
generally, we can conjecture that only the boundary conditions such that the variation of 
the function on the boundary can remain bounded will lead to the usual thermodynamic 
limit. This is essentially what is proven in [13] in the case A = 0. 

6. Conclusion 

In this work, we have computed explicitly the large N asymptotic behavior of a 
N X N determinant which plays the role of partition function of the six-vertex model 
with domain wall boundary conditions. This gives rise to particularly simple expressions 
for the bulk free energy of this model (Eqs. (4.10) and (4.18)). One important question 
is to physically interpret the discrepancy of the bulk free energy found when comparing 
domain wall and periodic boundary conditions of the six-vertex model, which is somewhat 
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contrary to standard lore on the thermodynamic hmit of statistical models. Some clues 
were given in the previous section, where various equivalences were discussed. In particular 
it was shown how "generic" fixed boundary conditions for the six- vertex model do not lead 
to a well-defined thermodynamic limit. It would be useful to make these arguments more 
rigorous. Also, it would be most interesting to find a more quantitative description of the 
non-translational invariance created by the boundary conditions, and in particular to prove 
a "generalized arctic circle theorem" for any value of the parameter A of the six-vertex 
model. 
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